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Combustion noise from turbofan engines has become important, as the noise from sources like the fan 
and jet are reduced. An aligned and un-aligned coherence technique has been developed to determine 
a threshold level for the coherence and thereby help to separate the coherent combustion noise source 
from other noise sources measured with far-held microphones. This method is compared with a statis- 
tics based coherence threshold estimation method. In addition, the un-aligned coherence procedure at 
the same time also reveals periodicities, spectral lines, and undamped sinusoids hidden by broadband 
turbofan engine noise. In calculating the coherence threshold using a statistical method, one may use 
either the number of independent records or a larger number corresponding to the number of overlapped 
records used to create the average. Using data from a turbofan engine and a simulation this paper shows 
that applying the Fisher z-transform to the un-aligned coherence can aid in making the proper selection 
of samples and produce a reasonable statistics based coherence threshold. Examples are presented 
showing that the underlying tonal and coherent broad band structure which is buried under random 
broadband noise and jet noise can be determined. The method also shows the possible presence of 
indirect combustion noise. [DOI: 10.1121/1.3546097] 
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Resolution bandwidth, Hz, B e = A f= \/T d 
= r/NP= 11.71875 Hz 

Number of points contained in deliberate 
time delay shift interval, D = rtj 
Confidence interval 
Frequency, Hz 

Auto power spectral density function 
defined for non-negative frequencies only 
(one sided) 

Cross power spectral density function 
defined for non-negative frequencies only 
(one sided) 

Number of degrees of freedom 

Number of frequencies, f c /Af=NP/2 (2048) 

Number of frequencies used in simulation 

Corrected rotor speed, rpm 

Independent stationary Gaussian random 

time series 

Number of disjoint (independent) segments 
used in spectra estimates, n d = B e T tota i = 234 
Number of segments/blocks with overlap 
Segment length, number of data points per 
segment (4096) 

Probability confidence interval, % 

Sample rate, samples/s (48 000) 

Time, s 

Total record length, s («20 s) 

Record length of data segment i, Ttotai /«</, 
= = 4096/48 000 = 0.08533 s 
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The coherent output power function 
w transform function, n^y^) = tanh Ty^) 
Signal x time history 
Noise time series 

Signal y time history 
z transform function, z(x) = tanh _1 (^x) 
Signal x 
Signal y 

Noise time series mixing parameter 
Frequency step, 1 IT d , Hz (11.718) 

Sampling interval, l/r= 1/48 000), s 
Coherence threshold using a probability, P h 
and n records, n d 

Coherence threshold using the number of 
records, n, and a 95% confidence interval 
Estimated magnitude squared coherence 
(MSC) function 

Estimated un-aligned MSC function of 
signals 

Simulated coherence target value 
Deliberate time delay shift interval, r,= D/r, s 


I. INTRODUCTION 

Identification of noise sources in the core of a turbofan 
engine is difficult since the core noise signal consists of broad- 
band noise, tonal noise, and hydrodynamic pressure perturba- 
tions. This paper discusses measurements made by external 
far-field microphones and pressure sensors in the combustor. 
All the spectra and cross-spectra are estimated using Welch’s 
non-parametric method which is based on averaging multiple 
windowed periodograms using overlapping time sequences. 1 
The method is discussed in textbooks by Canday, 2 Stearns and 
David, 3 and by Manolakis el al. 4 
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Using these spectra and cross-spectra, the magnitude 
squared coherence (MSC) between a combustor sensor and 
a far-held microphone is calculated to measure the similar- 
ity of the amplitude variations at particular frequencies. 
Interpreting the coherence measurement was made difficult 
due to the presence of an external extraneous low fre- 
quency noise source, as discussed by Miles, 5 that reduces 
the measured coherence and required the calculation of a 
coherence threshold. The coherence threshold would be the 
coherence of two independent signals measured in the 
same manner which would indicate the lowest coherence 
measurable for two independent signals for a particular 
system. If the measured coherence is below the threshold, 
it is attributed to random noise. This paper presents the 
method used to establish the coherence threshold used for 
the turbofan engine data. Then results are presented, which 
show that the measured coherence between the combustor 
and the signal at a far-held microphone while small is 
above the coherence threshold. 

A statistical test for independence is discussed by Car- 
ter, 6 ' 7 Halliday et al ., 8 and Brillinger. 9 It provides an equa- 
tion to calculate the upper 95% confidence limit for the case 
when the measured functions are independent. The statistical 
test is a confidence interval, which is function of the number 
of independent samples (record segments), n d . The confi- 
dence interval is based on a statistical distribution of the 
magnitude squared coherence (MSC) estimates, y 2 y , com- 
puted from signals produced by two Gaussian random proc- 
esses that have been divided into n d independent segments 
having spectral densities estimated by Welch’s averaged 
periodogram method without overlapping segments. The test 
does not take into account the procedure for calculating the 
coherence such as the amount of sample overlap and does 
not account non-Gaussian signals containing tones. The sta- 
tistical test results are of interest because they provide a con- 
fidence interval independent of frequency. 

This statistical test for independence produces a 
threshold value that does not depend on the measurement 
procedure used. However, selection of the signal process- 
ing calculation procedure is the responsibility of the 
researcher. Consequently, a coherence threshold value may 
depend on the selection of the non-parametric coherence 
signal processing calculation procedure such as the Black- 
man-Tukey method discussed by Blackman and Tukey 10 
and by Manolakis et al 4 or Welch’s method discussed by 
Welch, 1 and by Manolakis et al. 4 It is also a function of 
the selection of a spectral window function as discussed by 
Canday, 2 Stearns and David, 3 and Manolakis et al. 4 
Another variable is the selection of the amount of segment 
overlap. In addition, the coherence threshold might depend 
on the statistical distribution of the time series and the 
presence of tones also. 

In an earlier paper. Miles 11 discusses a procedure that 
uses Welch’s method to calculate the cross-spectra of two 
time histories recorded simultaneously with the modification 
of deliberately translating one of them with respect to the 
other by a time delay interval, T| , that exceeds the sample re- 
cord length, T d , of the averaging periodogram spectrum ana- 
lyzer used. This procedure creates two distinct independent 


time histories that will have the same sinusoidal content. 
The measured coherence at frequencies where tones are not 
present will have a threshold value of two independent sig- 
nals. In addition, it is calculated on the identical measure- 
ment system used in calculating the aligned signal 
coherence. Since the translation does not change the phase 
relationships in the time sequence, all tonal information will 
appear. The procedure at the same time as it calculates a co- 
herence threshold also reveals periodicities, spectral lines, 
and undamped sinusoids hidden or concealed by broadband 
turbofan engine noise. This calculated coherence threshold 
will be designated as |y 2 (/)|. Moreover, the two time series 

processed can be identical. Then one produces a coherence 
of the auto-spectrum that shows the coherence threshold of 
the system and tones because all broadband random informa- 
tion is removed. 

As a consequence, the results of the statistical test for in- 
dependence can now be directly compared with an estimated 
coherence value produced by a simple signal analysis calcu- 
lation procedure. A surrogate data analysis procedure to cal- 
culate only the threshold coherence was studied by Faes 
et al., 12 and it is discussed in Sec. II. These measurement 
procedures endow the threshold coherence with a reality and 
may help identify subtle noise sources. 

The total turbofan engine data record length studied 
herein is about 20 s and engine tones are persistent through- 
out this interval. The tones are of interest so that some of the 
irregularities in the auto-spectrum can be identified and tones 
below the spectrum can be observed. Tracking spectral lines 
in signals with additive noise as a function of time from 
discrete-time observations is a well-known problem some- 
times called in the time series literature the identification of 
sinusoids-in-noise problem or the hidden periodicity prob- 
lem. Currently, problem is generally treated using a modern 
approach such as parametric spectral-estimation. One of the 
first signal processing textbooks discussing the classical and 
more modern approaches is by Candy. 2 A more recent signal 
processing textbook discussing the classical and more mod- 
ern approaches is by Manolakis et al 4 A mathematical 
approach to classical and modern signal processing topics is 
presented in book by Byrne. 13 The topic of sinusoids-in- 
noise is mentioned by Byrne. 1 ' 1 Several approaches to solve 
the problem of sinusoids-in-noise are discussed in the book 
by Quinn and Hannan. 14 This identification of sinusoids-in- 
noise problem which has applications to radar, communica- 
tions, audio, power system control, and medical fields is 
solved using very different methods than the one identified 
herein. The un-aligned coherence method is not designed for 
tone tracking. 

As discussed by Miles, 11 the un-aligned coherence 
method is related to the time delay bias error. It is well- 
known that if two time series are misaligned by a time delay 
interval, iq, the estimated MSC for this case, y 2 y , is reduced 
by the square of 1 minus the ratio of the time delay, tq to the 
data segment record length, T d 
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This relationship has been used to explain time delay bias by 
Bendat and Piersol. 15 The importance of eliminating time 
delay bias and realigning time series by displacing one of 
them with respect to the other to maximize coherence has 
been mentioned in books by Jenkins and Watts, 16 Priestley, 17 
and Brillinger 9 and in articles by Carter. 7 ' 16 The concept 
used herein involves using a deliberate time delay that 
exceeds the sample record length, for example, letting 
Ti « 1 .57',/ to find the coherence threshold, hidden periodici- 
ties, concealed spectral lines, and undamped sinusoids in tur- 
bofan engine signals. None of the above references mention 
this deliberate unalignment procedure. The only mention of 
its use appears in paper by Gross et a/. 19 in the medical liter- 
ature where it is used to determine a coherence threshold 
without any mention of the origin of the procedure and for a 
study where no mechanical tones are present. 

In the previous paper by Miles, 11 Welch’s averaged 
periodogram calculation procedure using a 50% overlap, as 
advised by Carter et a/., 18 was used to calculate the turbofan 
engine coherence. The overlap reduces the variance and as a 
consequence the number of samples of the time history used 
to calculate the coherence was n s = 2 n d where n d is the num- 
ber of independent samples or records. The statistical proce- 
dure used to calculate the threshold as derived is based on 
the number of independent samples, n d . The threshold value 
decreases as n d increases. In the previous paper by Miles, 11 a 
simulation was done, which indicated that for Gaussian dis- 
tributed random variables, the statistical test for the coher- 
ence threshold could use twice the number of independent 
records for a 50% overlap in signal processing rather than 
the number of independent records. Just as the choice of 
window function in a non-parametric spectral-estimation, 
calculation is a matter of judgment, so is the choice of n s in 
calculating the statistical coherence threshold. The present 
paper more closely examines this premise. It indicates how 
applying the Fisher z-transform -0 of the un-aligned coher- 
ence can aid in making the proper selection of samples, i.e., 
n s and produce a reasonable statistics based coherence 
threshold. The present paper presents results of a computer 
simulation and a turbofan engine test program. The results 
show that for the turbofan engine data set using the number 
of independent records, n s = n d produces a higher threshold 
which is more reasonable. However, the use of n s =2n d 
which produces a lower threshold is line for the computer 
simulation results which used a Gaussian distribution of ran- 
dom numbers. Results are presented, which show that while 
the coherence is low it is still above the higher coherence 
threshold based on using n s = n d and a small but discernable 
contribution of the total noise is due to combustion noise. 

Internal noise sources (i.e., core noise) become the sig- 
nificant contributors to overall turbofan engine noise during 
takeoff or approach when the fan and jet noise are reduced 
because of the forward velocity effects. The core noise con- 
sists of noise associated with combustor, compressor, and 
turbine noise. This paper discusses an extension of the three 
signal coherence technique to study the compressor and tur- 
bine noises and an extension of the coherent output power 
technique to study the combustor noise. The three signal co- 
herence technique was developed by Chung et air 1 and 
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Chung 22 for flow noise rejection. The basic formulation for 
the coherent output power spectrum is presented by Bendat 
and Piersol. 15 These methods depend on the measurement of 
the MSC function, y 2 r The coherence values are low so that 
to interpret these measurements, it is necessary to know the 
coherence threshold. 

The data analyzed is from a Pratt & Whitney PW4098 
turbofan engine (East Hartford, CT). The test was conducted 
as a part of the NASA Engine Validation of Noise Reduction 
Concepts (EVNRC) Program in 2001. Two pressure trans- 
ducers are mounted in the combustor and four far-field 
microphones are used. The combustor pressure transducers 
failed slowly during the test. Consequently, a sensor valida- 
tion analysis was conducted based on the aligned and un- 
aligned coherence method and is reported by Miles. 23 The 
PW4098 combustor pressure sensor cross-spectra amplitude 
and phase plots resembled the cross-spectra used in the 
acoustic modal analysis of an AVCO Lycoming YF102 com- 
bustor (Stratford, CT), installed in a ducted test rig. This 
YF102 test program was conducted by Karchmerr 4 This 
suggested that the observed structure in the coherence was 
due to a modal pattern. Consequently, using a model of the 
physics of pressure waves propagating in annular ducts as 
discussed by Tyler and Sofrin 25 in a treatment of axial flow 
compressor noise, a procedure discussed by Miles 26 was 
developed that uses the available combustor pressure sensor 
data to obtain a descriptions of the combustion modes in the 
annular combustor. Results from this work were used to 
select the signals to be analyzed and to interpret the results 
in an attempt to identify combustion noise in the far-field 
data. Work was then done that showed that by comparing 
aligned coherent power spectra calculated using the esti- 
mated time delay with the un-aligned coherent power spectra 
calculated with no time delay one can identify combustion 
noise. In addition, it was shown that conventional statistical 
methods could be used to detect correlated combustion noise 
and create a one parameter measure of its presence or ab- 
sence. It was shown that the correlated combustion noise is 
only a small part of the total noise and that the plane wave 
combustion mode is apparent at low engine speeds. The 
aligned and un-aligned coherence method used in this work 
is discussed herein. Results are illustrated with examples 
from studies made of the Pratt and Whitney P&W 4098 tur- 
bofan engine data. 

The aligned and un-aligned coherence procedure of 
Miles 11 was used to develop a method to validate combustor 
coherence measurements by Miles 2 ’ and to find the coher- 
ence threshold in a study of turbofan engine combustor an- 
nular duct modes by Miles. 26 The method was also used to 
find the coherence threshold in analyzing core noise meas- 
urements made in a smaller Honeywell TECH977 turbofan 
engine (Phoenix, AZ). The smaller engine test results are 
reported in papers by Miles discussing a direct and indirect 
combustion noise source separation research program. 27 -9 
The core noise results discussed herein will be interpreted in 
terms of these more recent results. 

The deliberate decorrelation of signals was used in a 
study of flow-induced sound of wall-mounted cylinders with 
different geometries by Becker et al., 30 a study of source 
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separation from multiple microphone measurements in the 
far-held of a full scale aero engine by Mendoza et al } 1 and 
in an engine validation of noise and emission reduction tech- 
nology report by Weir. " 

For completeness, in Sec. II additional studies of coher- 
ence threshold calculation are discussed. Then the experi- 
ment and core noise related coherence definitions are 
discussed. Next the statistical estimate of the coherence 
threshold is discussed. Then the aligned and un-aligned co- 
herence method and results are presented. In presenting the 
results, first some simulation results are presented. Then, 
some experimental Pratt & Whitney PW4098 turbofan 
engine results are analyzed. Next, the results are discussed 
and conclusions are presented. 

II. OTHER THRESHOLD COHERENCE STUDIES 

A slightly different aligned and un-aligned coherence 
procedure was used by Gross et al. 19 in a study of cortico- 
muscular synchronization during isometric muscle contraction 
in humans as revealed by magnetoencephalography. After cal- 
culating the aligned and un-aligned coherence, the confidence 
limit yl n (\ — a.) is estimated as a value under which lies 
(1 — a) 100% of the values of the estimated MSC function in 
the frequency range used. Gross used a 50% overlap. 

The procedure used by Gross et al. 19 was examined by 
Bortel and Sovka 33 in a study of an approximation of the sta- 
tistical distribution of the MSC function estimated with seg- 
ment overlapping. They are unhappy with this estimation 
procedure. First, since it is based on the estimation of a 1 — a 
quantile, which is an imprecise operation for which there are 
many procedures, it can lead to different confidence limit 
estimates. Second, since various time delays D produce dif- 
ferent MSC functions this can also lead to different confi- 
dence limit estimates. Bortel and Sovka 13 go on to develop 
an approximation to the statistical distribution of the MSC 
function with segment overlapping and to develop a proce- 
dure for confidence interval estimation. Bortel and Sovka 11 
suggest that the value n d used in calculating the statistical 
coherence threshold corresponds to the number of degrees of 
freedom, L. They then speculate that the statistical distribu- 
tion of his time series is similar to that used in developing 
the statistical threshold and that details like the amount of 
overlap and the window can be accounted for by the proper 
choice of the number of degrees of freedom, L. Next, estima- 
tion equations are derived for L. 

A surrogate data analysis procedure studied by Faes 
et al} 1 has been used in biomedical engineering studies to 
develop a threshold for the coherence of two uncorrelated 
time histories. The paper by Faes et al. examines three pro- 
cedures for generating surrogate series. The method of surro- 
gate data generates a set of surrogate series mimicking some 
properties of the original series but being otherwise 
uncoupled by for example randomly shuffling sample 
blocks. Using a coherence of the surrogate data, the confi- 
dence limit y”„(l — a) is estimated as a value under which 
lies (1 — a) 100% of the values of the estimated MSC func- 
tion in the frequency range used. Consequently, the surro- 
gate data procedure also depends on the method used for 



FIG. 1. (Color online) Simulation result summary. 

calculating the coherence such as the amount of sample 
overlap. However, it cannot identify tones since the random 
shuffling of data sample blocks should remove periodic in- 
formation. In addition, the method has a similar lack of pre- 
cision as aligned and un-aligned coherence method, since 
different procedures for generating the surrogate series might 
give different results and many procedures are available for 
calculating the contidence limit. 

In this paper, the un-aligned and un-aligned coherence 
procedure and the threshold statistical procedure using L = n d 
and L = 2 n d are examined and compared. The surrogate data 
analysis procedure studied by Faes et al} 1 and the procedure 
suggested by Bortel and Sovka 11 were not examined. 

III. EXPERIMENT 

To demonstrate the usefulness of the proposed aligned 
and un-aligned coherence method, pressure measurements 
made in a Pratt & Whitney PW4098 combustor will be used. 
The measurements were made in a study of aircraft turbofan 
engine core noise conducted as part of the NASA EVNRC 
Program at test stand C-l 1 in West Palm Beach, Florida. The 
engine had an attached acoustic inflow control device. A pho- 
tograph of the test site appears in Fig. 1 of Ref. 11. The com- 
bustor was instrumented with two pressure sensors. One 
water cooled combustor pressure sensor was at 127° (sensor 
1) and the other was at 337° (sensor 2). A photograph of the 
mounted hardware is shown in Fig. 2 of Ref. 11. Combustor 
pressure sensor angles are measured clockwise from top dead 
center viewed from the rear or exhaust section of the engine. 
The spectral estimate parameters are shown in Table I. The 
signal processing procedure used is discussed in Sec. IV. The 
far- field microphones used herein are on a 150 foot radius at 
angles of 100° (sensor 3), 110° (sensor 4), 120° (sensor 5), 
and 130° (sensor 6) measured from the inlet. The acoustic 
arena and microphone array are shown in Fig. 3 of Ref. 11. 

IV. CORE NOISE STUDY SYSTEM CONCEPTS 

The procedure discussed herein involves the following 
items: 

(1) The estimated MSC function, y 2 xy [f)'. 

First, the MSC function will be discussed. Then, the 

estimation procedure will be discussed. 
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— Simulated aligned coherence 
— O- Noise unaligned coherence 

Lower confidence interval 

upper confidence inverval 

Target cohenence, 0.03 

Noi S e coherence threshold, 0.006 (n s = 500) 

Noise coherence threshold, 0.012 (n<j = 250) 


TABLE I. Spectral estimate parameters. 


Parameter 

Value 

Segment/block length, i.e.. Data points 

4096 

per segment/block, NP 


Sample rate, r, samples/s 

48 000 

Segment length, T d — NP/r, s 

0.08533 

Sampling interval, At =l/r, s 

2.0833 x 10~ 5 

Frequency step, Af— B e = \/T d, Hz 

11.718 

Upper frequency limit, f c = 1/2 At = r/2, Hz 

24000 

number of frequencies, Ly —fJAf— NP/2 

2048 

Time delay to far field, % = 6323/48000, s 

0.1317 

number of independent segments/blocks 

234 

Overlap 

0.50 

Sample length, s 

20 

Spectral window 

Rectangular 



Norm quantiles 



w(Coherence) 


FIG. 2. (Color online) Plots show target coherence, measured aligned coher- 
ence, the confidence interval, the un-aligned noise coherence, the estimated 
upper 95% confidence limit threshold for a 50% overlap y% nyn (f, Pi = 0.95, 
n s = 500) ~ 0.006, and based on the assumption of independent records the 
estimated upper 95% confidence limit threshold y^ y ^{f : Pj = 0.95, 
rid = 250) ^ 0.012. (a) Simulation using target coherence {y^ nyn \ T = 0.03. (b) 
Norm quantile comparison of un-aligned noise coherence, (c) Histogram of un- 
aligned noise coherence. 


fxy(f) = 


\G xy (f )\ 2 

G xx [f)G yy (f) 


( 2 ) 


The MSC function is discussed by Bendat and Pier- 
sol , 15 Carter , 7 Halliday et al. , 8 and Jenkins and Watts 16 
where it is called the squared coherency spectrum. 

The MSC function estimation procedure used herein 
are based on signal processing algorithms developed 
by Stearns and David.’ They refer to their spectral- 
estimation code as using a periodogram averaging method. 
The algorithm is also referred to a being based on Welch’s 
method (Welch 1 ) and as a weighted overlapped segment 
averaging (WOSA) method (Carter et a/. 18 ). 

Two additional characteristics of the MSC estimation 
process is its sensitivity to the alignment or translation of 
the time series and the inability to observe tones mixed 
with broadband noise in the spectrum. These points will 
be discussed later. One of the objectives of this report is 
to point out how deliberate misalignment of the time se- 
ries can enable one to observe tones in the spectrum. 



-A- -Aligned coherence 
--o- - Unaligned coherence 

Threshold coherence, 0.006 

Threshold coherence, 0.012 

Upper confidence interval 

Lower confidence interval 


The MSC function, y 2 xy (f), calculated using the com- 
plex cross-spectral density, G xy (f), and two auto-spectral 
densities at frequency/, G xx (f) and G yy (f), is defined by 


FIG. 3. (Color online) Aligned and un-aligned coherence between combus- 
tor pressure sensors 1 and 2. The broken horizontal line is (/, 
P[ — 0.95, rid = 250) ^ 0.012. The unbroken horizontal line is yx n y n (fi 
P I = 0.95, n s = 500) « 0.012. N1 Corr= 1622. 
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(2) The spectrum based on the three signal magnitude coher- 
ence function developed by Chung et al. 21 and Chung 22 for 
flow noise rejection is (also see Appendix B of Ref. 1 1 ) 

r (f\ — r \ m% m v)\ ... 

a (fl=G ™ (7) ' <3> 


(3) The coherent output power spectrum function discussed 
by Bendat and Piersol 15 (also see Appendix C of Ref. 1 1) 

Ws F s F (f) = f xy (f)G yy (f). (4) 

V. TRADITIONAL VIEWPOINT 

First, the traditional view on measuring coherence will 
be presented to pinpoint the “blind spots.” A good place to 
start is a statement about the proper analysis procedures to 
obtain a coherence measurement as stated by Piersol: 14 

(1) The instrumentation is calibrated to eliminate relative 
phase errors in the measured signals. 

(2) The received signals are aligned as required to eliminate 
time delay bias errors in the coherence calculation. 

(3) The number of averages n d is sufficiently large to suppress 
small sample bias errors in the coherence calculations. 

VI. STATISTICAL ESTIMATES OF COHERENCE 
A. Confidence interval 

A significant limitation to the application of the MSC 
estimate is that for coherence values smaller than some fairly 
large value (like 0.1) and for large sample sizes, the confi- 
dence interval is not available. For values greater than some 
small value (like 0.1), the confidence interval is only given 
approximately (Bendat and Piersol. 2 Brillinger, 9 and Car- 
ter 7 ' 18 ). A procedure calculating the exact confidence interval 
is discussed by Wang et al , 36 ' 37 and by Miles. 11 Conse- 
quently, if one calculates the upper and lower confidence 
interval one can determine the significance of the estimation 
and the significance of the differences of MSC values at dif- 
ferent frequencies. Confidence limits are shown in Fig. 4 of 
Ref. 1 1 for n d = 600. The confidence limits are shown on lin- 
ear plot as they are usually shown in Fig. 4(a) of Ref. 11. The 
confidence limits are plotted in a more realistic manner in 
Fig. 4(b) of Ref. 11 on a logarithmic plot, which demon- 
strates clearly that for a given value of n d , the confidence 
interval becomes larger as the value of y 2 xy is reduced. 



norm quantiles 

FIG. 4. (Color online) Quantile comparison plot of unaligned coherence 
between combustor pressure sensor (1) and (2) with normal distribution. 
The broken lines give a pointwise 95% confidence envelope around the 
fitted solid line. The upper horizontal line at w(\/0.012) = tanh -1 (0.10954) 
« 0.109 represents a statistical coherence threshold of y\ yn (f, Pi 
= 0.95, nd = 250) « 0.012. The lower horizontal line at w(\/0.006) 
= tanh -1 (0.07746) ^0.0776 represents statistical coherence threshold of 
f Xnyn (f, P, = 0.95, n s = 500) w 0.006. N1 Corr= 1622. 


certain probability. The threshold value of this confidence 
interval will be called E\. 

The conditional probability of the estimated value of the 
MSC, |y| 2 , given the true MSC is zero, \y\ 2 = 0, and the number 
of independent segments/blocks, n d , can be calculated from as 


p{\y\ 2 \nd, |y | 2 = o) = (n d - i)(i - lyl 2 )" 7 2 . (5) 

This equation is given by Brillinger 9 and is discussed by 
Miles. 11 

Then the probability, P h that v 2 is in the interval 
0 < |y | 2 < Ei can be calculated as 


C = 
P,= 


p{\f\n d , |y | 2 = o)d\y\ 2 


(, n d - 1)(1 -C)" d 2 dC 


= -(i -c) 


id - 1 


= -(1 -E,) nd ~ l + 1. 


( 6 ) 


B. Statistical estimate of coherence threshold of 
independent processes 

In considering the MSC of two stationary Gaussian ran- 
dom processes, it is often stated that the MSC will be zero 
when the processes are independent and do not have a linear 
time-invariant relationship. However, the estimated MSC of 
two independent processes is never zero but has some small 
numerical value at each frequency. The only statement that 
can be made is that if the true MSC value is 0, |y| = 0, the 
expected value is in a certain confidence interval with a 


Thus, the cumulative distribution derived using the probabil- 
ity density function of the coherence derived by Goodman 18 
has been used to determine a By percent confidence interval, 
if the true MSC is zero. This is given by 

y 2 Xnyn if,Pi,n d ) = E, = 1 - (1 - p,) l P nd ~ l \ (7) 

where n d is the number of independent segments used in the 
spectrum calculations and the interval is between 0 and 
y 2 ( f,P[,n d ). This relationship is tabulated in Table III by 
Miles 11 for probability values, P h of 80%, 90%, 95%, and 
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99% and is shown in Fig. 1. This statistical test for independ- 
ence using Eq. 7 is discussed by Carter, 6 ' 7 who proposes using 
y 2 (/, P/. nd) as a threshold value below which coherence 
estimates are not presented to a human decision maker look- 
ing at an extensive volume of data collected for brain wave 
analysis or sonar or radar signal analysis. This equation is also 
given by Gosselin, 39 Halliday et al., s and Brillinger. 9 Halliday 
suggests one should plot the value of y 2 (f . Pi, no) calculated 
from Eq. (7), i.e., under the hypothesis of independence on 
each coherence plot. He states that estimated values of coher- 
ence lying below this line can be taken as evidence to indicate 
that the coherence is zero at that frequency. This opinion is 
shared by Nuttall 40 who thinks this is the standard result for a 
true coherence of zero. However, the author thinks it may 
indicate that not enough segments/blocks have been used to 
determine the coherence. One can state with more confidence 
that values above the y 2 (f. Pi, nd) line are not from two in- 

dependent stationary random Gaussian series. Note that this 
statement is valid independent of the magnitude of 
y 2 (/, P[, nd). In this paper, we examine whether one should 
calculate the coherence threshold using 

fx n ytf’ p i’ n s) = 1 - (1 -^/) 1/( " I_1) , 


where with a sample overlap of 50% n s = 2 n d . In the figure 
legends, the notation for the threshold coherence drops the 
accents and becomes simply y 2 n («) where the n value in 
parentheses indicates one is using the number of independent 
samples n d or the number of overlapped samples 2 n d and 
the confidence interval probability, P t , is 95%. 


VII. COHERENCE TRANSFORM 

When the transformation 


z = l i n A 1 + yg ) N 

= tanfT 1 (v/A") 


(8) 


(termed the Fisher z-transformation 20 ) is applied to the ordi- 
nary correlation coefficient, one obtains a variable which is 
normally distributed to a very high degree of accuracy. Con- 
sequently, the random variable w 



/(l + 7xy)\ 

l(i -y»)J 


= tanh *.( iy 


(9) 


is to a close approximation, Gaussian with a distribution 
N(/.i, cr 2 ). 41 The wiy^) transforms of the estimated nonaligned 
coherence functions calculated herein will be compared to 
normal distributions using histograms and quantile-quantile 
(Q-Q) plots, which are the standard statistical analysis tools. 42 


VIII. COHERENCE AND SIMULATION METHODS 
A. Simulation procedure 

Many assessments of the validity of the statistical esti- 
mates of coherence have been made using simulations. The 


probability density function of the coherence derived by 
Goodman 8 was derived for Gaussian data. Foster and 
Guinzy 43 tested this distribution by means of Monte Carlo 
experiments for validity and robustness (insensitivity to the 
Gaussian assumption) and it passed the tests. Coherence 
function bias and confidence intervals were studied using 
Monte Carlo methods by Benignus. 44 

The probability density function of the coherence 
derived by Goodman 38 does not take into account the use of 
overlapping data. Empirical simulation results by Carter 7 
and Carter et al . 1 8 show a decrease in bias and variance of 
the estimator with increasing overlap and suggest a 50% 
overlap as being highly desirable. 

Furthermore, since the probability density function of 
the coherence derived by Goodman 38 is independent of the 
computer code time series analysis implementation compar- 
ing bias and variance simulation results can provide some 
confidence in the analysis procedure. Simulations were done 
for a range of target coherence values and a range of inde- 
pendent segments/blocks using 


x„(t) = « 3 + am, 
y n (t) = « 3 + om 2 , 

z„(0 = «3, (10) 

where « 1; n 2 , and « 3 are independent Gaussian random time 
series. Results will be presented in Figs. 1 and 2. 

The time series are selected so that 


Gn. i'i ; G n | R] c„ 2 „ 2 G ss , 

1 ‘ n 2 G ) ; H 1 if: : fi ; 0. 

Taking the Fourier transform of x„ yields 

Xx»(f)=Xn 3 (f) + 0tX ni (f). (11) 

The spectrum G XnXn is then 


=x Xn (f)x; n (f) 

= (X n ,(f) + xX ni (f))(X n3 (f) + aXimmy 
T a Gn,n, • 

We then also have G yn y n = G„ 3 „ 3 + ocG n2 „ 2 and 


G XnX „ 

— 1 

G x „z n 

— C"« 3 A ?3 i 

Gy„z n 

— ^rt 3«3 1 

Gz n z n 

— ^/Z 3«3 • 


Then the target MSC is a function of a given by 




G W ,G W , (1 + a 2 ) 2 


( 12 ) 


and 


{yL„(f)} T = {f yA (f)} T = 


(1 + a 2 ) ’ 


(13) 
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These expressions are independent of frequency. A similar 
approach is used to create target simulations by Carter 
et a/. 18 and Guinzy. 43 

B. Simulation results 

Results of a typical simulation using a target coherence 
of { y 3 (/)} = 0.03 are shown in Fig. 2(a). The simulations 
are done with n d = 250 with a 50% overlap, so n s = 500. The 
coherence is low, so a logarithmic coherence scale is used in 
Fig. 2(a). The target coherence is shown by the heavy dashed 
line. The plot shows the estimated or measured coherence, 
y 3 y (f) ’ as triangles and an estimated un-aligned signal coher- 
ence, |y 3 v (f) | as circles for n s = 500. The upper and lower 
confidence interval is shown as dashed lines. In addition, the 
lower horizontal statistical coherence threshold line 
yl n y n (f, Pi = 0.95, n s = 500) ~ 0.006 and the higher hori- 
zontal statistical coherence threshold line y 3 (/'. P t = 0.95, 
n s = 250) ~ 0.012 are shown. For this case, the estimated 
un-aligned signal coherence, |y 3 (f)\, is less than the lower 

statistical coherence threshold (y 3 (f 1 P = 0.95 ,n s = 500) 
ss 0.006) and the measured coherence is well above both sta- 
tistical coherence thresholds. Consequently, the coherence 
while low is significant. The lower statistical coherence 
threshold line traces the tops of the estimated un-aligned co- 
herence curve. The upper statistical coherence threshold line 
at y 3 y (f,Pi = 0.95, n s = 250) « 0.012 is at a value greater 
than that of the estimated un-aligned signal coherence, 
|y 3 (f ) |, which has values less than 0.006. 

The distribution of the un-aligned noise coherence plot- 
ted in Fig. 2(a) was examined to determine if it had a normal 
distribution using two procedures. A norm quantile compari- 
son plot of the Fisher z-transformed un-aligned noise coher- 
ence using Eq. (9) is shown in Fig. 2(b). The broken lines 
give a pointwise 95% confidence envelope around the fitted 
solid line. Horizontal lines corresponding to Fisher z-trans- 
formed un-aligned threshold coherence values of 0.006 and 
0.012 are also shown. The upper horizontal line in Fig. 2(b) 
at w(V 0.012) = tanlT 1 (0. 1095) « 0.109 represents an esti- 
mate of the upper 95% threshold coherence confidence limit 
of the un-aligned coherence based on 250 independent 
samples. The lower horizontal line at w(a/0.006) 
= tanh -1 (0.07746) « 0.077 represents an estimate of the 
upper 95% threshold coherence confidence limit of the un- 
aligned coherence based on 500 samples. This line is seen to 
provide a better coherence threshold for the 95% confidence 
interval and to be in good agreement with the estimated un- 
aligned signal coherence values, |y 3 v (/)|, shown in 
Fig. 2(b). The coherence threshold calculated with n d =250 
is larger than necessary. 

A corresponding histogram plot of the Fisher z-trans- 
formed estimated un-aligned signal coherence values, 
|y 3 (/') is shown in Fig. 2(c). Due to a lack of observations 

only four bins are used. Also shown are non-parametric 
kernel-density estimates for the distribution of coherence, 
using the default bandwidth (heavier line), which produces a 
smoother density estimate, and half the default bandwidth 
(lighter line), which produces a rougher density estimate. A 
normal distribution is shown as a dashed curve. A one- 


dimensional scatter plot is shown at the bottom of the graph, 
using a vertical bar as the plotting symbol, with the horizon- 
tal location given by coherence. On this plot vertical lines 
indicate estimates of the upper 95% statistical threshold co- 
herence confidence limit calculated with n d = 250 which pro- 
duces a line at h’( 0.012)« 0.109 and «, = 500 which 
produces a line at w(0.006) ~ 0.077. The line produced with 
n s = 500 is again seen to produce a better statistical coher- 
ence threshold value. 

The results of a great many simulations are shown in 
Fig. 1 where average MSC values, y 3 V , are presented as a 
function of the number of independent segments/blocks, n d . 

Each simulation resembles the one shown in Fig. 2(a) 
and is designed for a particular target value. We use as a 
measure of convergence to a particular target value 


1 


n— 1 


y% - „ ^ 

i=0 


fxyifi)' 


(14) 


Simulations with a range of target MSC values calcu- 
lated using Eqs. (9) and (10) are shown. The symbols trace 
each target value. The target values are given in the legend. 
In addition, simulation results for the mean value when two 
independent simulated random Gaussian time histories are 
used to calculate y 3 „ (f) are shown as black triangles and 

A "l/fi 2 

identified as y^, = 0 in the legend. 

The results show that a coherent pair of signals will have 
a coherence that resembles the coherence of random noise if 
the number of segments/blocks, n d , is too small. For example, 
for n d = 100, one would see an indication of MSC = 0.01 but 
not of MSC = 0.0025. Consequently, the value n d controls not 
only the bias and variance of the coherence but also its exis- 
tence. Consequently, for a given measurement of coherence 
using a given value of n d one can say that a coherence may not 
be observed if there is a P percent probability that it is less 
than yl nyn (f,P h n d )=E, = 1 - (1 — p I ) 1 l^ td ~ v > Eq. (7). This 
curve is identified as y = 1 - 0.05* in the legend box. This 
curve is shown in the figure as a thin line for a 95% confidence 
interval. In addition, the y 3 v = l/n d line is shown as a diagonal 
black line identified as y = x in the legend. This line corre- 
sponds to the expected or mean value of |y|~, for |y|^ v = 0 as 
discussed in Ref. 10 [see Eq. (12)]. This curve also corre- 
sponds to the one Tukey 45 suggested as discussed in Ref. 11. 
The estimated mean value of |y|^ (f) shown by the black tri- 

angles is close to this line. However, the 95% confidence inter- 
val provides a better number to be larger than, when one 
wishes to identify a small coherence value as being significant. 

The simulations showed that with a Gaussian distribu- 
tion, the coherence threshold calculated using the un-aligned 
coherence with a 50% overlap is similar to the statistical co- 
herence threshold calculated with n s = 2 n d . 


IX. ALIGNED AND UN-ALIGNED COHERENCE 

The method of aligned and un-aligned coherence is 
based on the following two characteristics of the average 
periodogram process when it is used to calculate the MSC: 

(1) When the deliberate time delay shift interval, t u exceeds 
the sample record length, T d , the coherence of the 
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random process, y 2 xy (f), is the coherence of two inde- 
pendent random signals |y 2 (f) |. 

(2) If the signals contain tones using a deliberate time delay 
shift, t | , greater than the sample record length, T d , does 
not change the coherence of the tones and permit the 
estimation of hidden periodicities, concealed spectral 
lines, and undamped sinusoids. These items are gener- 
ally referred to as tones herein. 

The deliberate time delay shift interval is sometimes 
expressed herein in terms of the number of data points con- 
tained in the time delay shift interval, D, which is the product 
of the sample rate, r, and the time delay shift, ti, ( D = r Ti). 
Spectral estimate parameters used herein are shown in Table I. 
For the data processed herein, the length of the sample record 
or segment interval, T s , is T d = NP/r = 0.08533 s since the 
number of data points contained in a segment/block is 
NP = 4096 and the sampling rate, r, is r = 48 000 samples/s. 
The necessary condition to get nonalignment can be expressed 
as deliberately displacing one time series with respect to 
another by shifting it an amount D such that D > NP or by dis- 
placing it in time by an amount Ti where t\ > T d . 

X. EXPERIMENTAL RESULTS 

A. Combustor pressure sensor coherence 

Logarithmic plots of aligned coherence, y 2 (f) (shown 
by filled triangles), and un-aligned coherence, \y 2 y (f)\ 
(shown by circles), between the two combustor pressure sen- 
sors are shown in Fig. 3 for a typical test conditions (N1 
Corr= 1622). Also shown are the confidence intervals for the 
aligned coherence. Note that as indicated in Fig. 4(b) of Ref. 
1 1 , the confidence interval is narrow for coherence values 
near 1 and becomes much larger for lower coherence values. 
In addition, the lower horizontal statistical coherence thresh- 
old line f Xny Jn s ) = y\ ny (f.P, = 0.95, n s = 500) « 0.006 
and the higher horizontal statistical coherence threshold line 
y 2 n y n { n d) = y 2 n y n (f,Pi = 0.95, n s = 250) « 0.012 calculated 
from Eq. (7) are shown. 

As previously mentioned, the PW4098 combustor pressure 
sensor cross-spectra amplitude and phase plots resembled the 
cross-spectra used in the acoustic modal analysis of an YF102 
combustor installed in a ducted test rig, which was conducted 
by Karchmer. 24 This suggested that the observed structure in 
the coherence was due to a modal pattern. Consequently, using 
a model of the physics of pressure waves propagating in annu- 
lar ducts as discussed by Tyler and Sofrin 25 in a treatment of 
axial flow compressor noise, a procedure was developed that 
uses the available combustor pressure sensor data to obtain a 
descriptions of the combustion modes in the annular combus- 
tor. The model is discussed by Miles. 26 Results from this work 
were used to select the signals to be analyzed and to interpret 
the results in an attempt to identify combustion noise in the far- 
field data. These broad peaks and dips in the coherence due to 
combustion duct modes are seen in Fig. 3. The coherent plane 
wave mode as determined by a model analysis is apparent in 
Fig. 3 between 0 and 200 Hz at 1622 rpm (N1 Corr.). 

Note in Fig. 3, the aligned coherence, y xy (f) (shown by 
diamonds) and the un-aligned coherence values \y 2 y (f)\ 


(shown by circles) are well separated. The aligned coherence, 
y 2 y (/), is well above both of the statistical threshold coherence 
horizontal lines [y 2 (f) and y 2 (fid)] and above the esti- 

mated un-aligned coherence, y 2 y (f) |, except in a region 
between 400 and 500 Hz. This bump appears to be a tone 
since it appears in both the aligned and un-aligned coherence. 

The same two procedures that were applied to the simu- 
lated time series are applied to the engine noise time series anal- 
ysis to determine if the un-aligned noise coherence shown in 
Fig. 3 has a normal distribution. A norm quantile comparison 
plot of the Fisher z-transformed un-aligned noise coherence 
using Eq. (9) is shown in Fig. 4 for a typical engine operating 
condition (N1 Corr = 1622 rpm). The broken lines give a point- 
wise 95% confidence envelope around the fitted solid line. Hori- 
zontal lines corresponding to Fisher z-transformed un-aligned 
noise coherence threshold values of 0.006 and 0.012 are also 
shown. The upper horizontal line in Fig. 4 at h’v/ 0.012 
= tanh '(0.1095) « 0.109 represents an estimate of the upper 
95% confidence limit of the un-aligned coherence threshold 
based on 250 independent samples. The lower horizontal line at 
w V0. 006 = tanh -1 (0.07746) « 0.077 in these figures repre- 
sents an estimate of the upper 95% confidence limit of the un- 
aligned coherence threshold based on 500 samples. However, 
for this turbofan engine data case the upper horizontal line based 
on n d = 250 provides a better coherence threshold. 

The corresponding histogram plot of the Fisher z-trans- 
formed un-aligned noise coherence is shown in Fig. 5 for the 
N1 Corr= 1622 rpm case. Also shown are non-parametric 
kernel-density estimates for the distribution of coherence, 
using the default bandwidth (heavier line), which produces a 
smoother density estimate, and half the default bandwidth 



FIG. 5. (Color online) Histogram of un-aligned coherence between combus- 
tor pressure sensors 1 and 2. Vertical lines w(\/0.006) = tanh -1 (0.07746) 
= 0.0776 and w(\/0.012) = tanh -1 (0.10954) = 0.109 represent statistical 
coherence thresholds of y\ n y n if, Pi = 0.95, n s = 500) « 0.006 and y^ n y n 
(/, Pi — 0.95, rid = 250) ~ 0.012 Non-parametric kernel-density estimates 
for the distribution of coherence, using the default bandwidth (heavier line) 
and half the default bandwidth (lighter line). A one-dimensional scatter plot 
at the bottom of the graph uses a vertical bar to indicate the horizontal posi- 
tion of each coherence magnitude value. N1 Corr= 1622. 
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FIG. 6. (Color online) Scatter plot of culled un-aligned coherence points 
from coherence between combustor pressure sensors 1 and 2. The coherence 
threshold based on n d = 250 is shown by the broken horizontal line at 
y 2 X]t y(f -Pi = 0.95, n ( i = 250) s» 0.012. The coherence threshold based on n, 
= 500 is shown by the unbroken horizontal line at X v (/ . / J / 0.95. 

n d = 500) ss 0.006. 


(lighter line), which produces a rougher density estimate. A 
normal distribution is shown as a dashed curve. A one- 
dimensional scatter plot is shown at the bottom of each 
graph, using a vertical bar as the plotting symbol, with the 
horizontal location given by coherence. On this plot, vertical 
lines indicate statistical estimates of the upper 95% threshold 
coherence confidence limit of the un-aligned coherence. 
Using n d = 250 produces a vertical line at w(0.012) « 0.109 
since y 2 (250) = 0.012 while using n s = 500 produces a 
vertical line at w(0.006) « 0.077 since y 2 y (500) = 0.006. 
Again, for this turbofan engine data case, the vertical based 
on n d = 250 provides a better coherence threshold. 

To investigate the frequency dependence of the engine 
coherence data with coherence greater than a threshold co- 
herence value of y 2 (n s = 500) = 0.006, these points were 
culled from a number of un-aligned coherence plots for a 
range of operating conditions and are shown in a scatter plot 
in Fig. 6. In addition, the lower horizontal statistical coher- 
ence threshold line 7^ >Vi (n s ) = 7x n yjf- = 0.95, n s = 500) 
« 0.006 and the higher horizontal statistical coherence 
threshold line y 2 ^ (n d ) = y 2 ^ (f, Pi = 0.95, n s = 250) 
~ 0.012 calculated from Eq. (7) are shown. This plot clearly 
indicates sinusoids-in-noise between 400 and 500 Hz. In 
addition, the broken horizontal line with a coherence thresh- 
old value y 2 ( n d = 250) = 0.012 is clearly a better choice 
of a coherence threshold. With this decision, most of the 
scattered culled points are interpreted as being due to statisti- 
cal random fluctuations and not tones. 

Comparing the engine test un-aligned results shown in 
Fig. 3-5 with the simulated noise un-aligned results shown in 
Figs. 2, it is obvious that using the value of the number of dis- 
joint (independent) segments, n d , provides a better estimate 
of the coherence threshold for the engine data. However, 
using the number of segments with overlap, n s , provides a 
good threshold estimate for the data simulated using a Gaus- 
sian distribution of values. In the next few sections, turbofan 
engine core noise measurements using the three signal coher- 
ence technique and the coherent output power calculation are 
examined using these statistical coherence threshold values. 



-- G 3 3=|Gi3(D=0)||G23(D=0)|/|G 12 (D=6323)| unaligned 
— G 3 3=|G 1 3(D=6839)||G23(D=6134)|/|G 12 (D=0)|) aligned 
' G 33 Ynn( n d = 250) 

G 33 Ynn( n d = 500) 


FIG. 7. (Color online) Three signal coherence technique using one far-held 
microphone and two combustor pressure sensors. The dashed line is based 
on an estimate of the upper 95% confidence limit of the un-aligned coher- 
ence threshold using n d = 250 (y^(f, Pi = 0.95, n d = 250) sa 0.012). The 
heavy line is based on an estimate of the upper 95% confidence limit of the 
un-aligned coherence threshold using n s = 500 (y) (/. P, = 0.95, 
n s = 500) s« 0.006). 


B. Three signal coherence technique 

The three signal coherence technique 21 ' 22 is applied to 
noise from a P&W4098 turbofan engine using the aligned and 
un-aligned coherence functions. Typical results of applying 
the three signal and un-aligned coherence technique are 
shown in Figs. 7 and 8. Results calculated using signals from 
one far-field microphone and two combustor pressure sensors 
are shown in Fig. 7. Results calculated using signals from 
three far-field microphones are shown in Fig. 8. 

1. Three signal coherence calculated using one 
far-field microphone and two combustor 
pressure sensors 

The example shown is for 1622 rpm (N1 Corr.). The 
three signal aligned and un-aligned coherence calculated 
using the microphone at 100° (sensor 3) and the two com- 
bustor pressure sensors (sensors 1 and 2) are used to show 
the presence of combustion noise in the region from 0 to 200 
Hz in Fig. 7. The results are shown so that they can be com- 
pared to auto-spectrum, G33, shown as the top curve com- 
posed of squares in Fig. 7. The aligned three signal 
spectrum, G(f) V3V} I aligned ’ shown by black circles is given by 


G (/% 3V3 


_r(f\ y(f)i3y(f)23 

aligned ~ )y 3 y 3 y(f) 


aligned • 


The un-aligned three signal spectrum, G(/) V3V J unaUgned , 
shown by triangles is given by 


G(f) v 


unaligned 


= G(f) 


J3B 


y(f)ny(f)23 

7(f) 12 


unaligned • 
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FIG. 8. (Color online) Three signal coherence technique using three far- 
held microphones. The dashed line is based on an estimate of the upper 95% 
confidence limit of the un-aligned coherence threshold using nj — 250 
(y^ n y n if i Pi = 0.95, rid = 250) ~ 0.012). The heavy line is based on an esti- 
mate of the upper 95% confidence limit of the un-aligned coherence thresh- 
old using n s — 500 (y% n y n (f, Pi = 0.95, n s = 500) » 0.006). 


The upper and lower coherence threshold spectra are also 
shown. The upper threshold G^y 2 : n (nd) is shown by the 
dashed curve. The lower threshold G^y 2 (n s ) is shown by the 
black curve. The un-aligned three signal spectrum 
G(fL J unaligned is below the coherence threshold at most fre- 
quencies except for the presence of tones. The aligned three 
signal spectrum G(f) \ aligned is above the coherence 
threshold curves in the region 0-200 Hz. This shows the 
presence of combustion noise. Results from a study of noise 
from a smaller Honeywell TECH977 dual-spool turbofan 
engine by Miles 27-29 suggest this noise is due to an indirect 
combustion noise source. 

The translation procedure will be discussed next. The 
number of translation points, D , used for each cross-spectrum 
so that the signals of interest would be in proper alignment is 
given in the legend. To calculate the aligned G 33 spectrum, a 
delay of 6839 time steps was used in calculating the aligned 
cross-spectrum G 13 to remove the time delay caused by the 
signal traveling between and microphone 3 and combustor 
sensor 1 . A delay of D = 6134 time steps was used in calculat- 
ing the G 23 aligned cross-spectrum. The delay times were 
selected by maximizing the coherence and making the plot of 
the cross spectrum phase angle as a function of frequency flat. 
The G 12 cross-spectrum was used as measured. To calculate 
the un-aligned G 33 spectrum a delay of 6323 time steps was 
used in calculating the cross spectrum G 12 , and the G 13 and 
G 23 cross-spectrum were used as measured. The D = 6839 and 
£> = 6134 time steps are in the range of the travel time from 
the turbofan to the microphones used. However, some attempt 
was made to find time intervals that maximized the coherence 
in the 0 to 200 Hz frequency band and those are used here. 

Tone noise appears in both the aligned and un-aligned 
poser spectrum. The tone at one BPF (first blade passing fre- 


quency) is especially noticeable. However, it is best shown 
in the un-aligned power spectrum where the broad band 
noise does not appear. Note that one can view tones approxi- 
mately 10 dB below the measured auto-spectrum in the three 
signal coherent aligned and un-aligned power spectrum. In 
the un-aligned three signal spectrum many tones previously 
unseen are now exposed to sight since the coherent broad- 
band noise has been removed. The aligned and un-aligned 
procedure solves the sinusoids-in-noise and hidden periodic- 
ity problem for turbofan engine noise. 


2. Three signal coherence calculated using far-field 
microphones 

Three signal aligned and un-aligned power spectra cal- 
culated using far-field microphones at 100° (sensor 3), 110° 
(sensor 4), and 120° (sensor 5) are shown in Fig. 8 . As with 
Fig. 7, the results are shown so that they can be compared to 
auto-spectrum, G 33 , shown as the top curve composed of 
squares in Fig. 8 . The plot clearly shows the tonal compo- 
nent due to the core and fan in the far-field noise spectrum, 
Gu 3 u 3 I unaligned, created using the un-aligned coherence (light 
diamonds) where 


G (/)m 3 h 3 unaligned ~G(f) 


y(/')3 4 y(/’)35 
7(f) 45 


unaligned • 


Above this spectrum is a spectrum G„ 3 „ 3 1 aligned created using 
the aligned coherence (dark dots) where 


G (/)u 3 H3 aligned — G(f) 


7(f)347(f)3S 
7(f) AS 


aligned • 


The upper and lower coherence threshold spectra are also 
shown. The upper threshold G 337 2 (n,j) is shown by the 

dashed curve. The lower threshold G^y 2 (n s ) is shown by 
the black curve. The un-aligned three signal spectmm 
g(/)„ 3 J unaligned is below the coherence threshold at most fre- 
quencies except for the presence of tones. The aligned three 
signal spectrum G(f) v v | aligned is far above the coherence 
threshold curves in the region 0-500 Hz. This shows the pres- 
ence of a strong low frequency coherent noise source. 

The aligned coherence values were used as measured. 
The un-aligned coherence values were found by delaying 
one signal 6323 time units with respect to the other. The jet 
noise is missing in the three signal aligned and un-aligned 
power spectra to the extent that the microphones are far 
enough apart that the jet noise signals they receive are inde- 
pendent in this frequency range. This feature of jet noise is 
exploited by the three signal coherence technique using three 
far-field microphones. The aligned and un-aligned spectra 
from 0 to 500 Hz are well separated showing the presence of 
a coherent noise source. 

However, this noise source could not be attributed to 
core noise. Furthermore, interpreting the coherence mea- 
surement used to obtain the coherent output power was 
made difficult due to the presence of an external extrane- 
ous low frequency noise source as discussed by Miles 5 
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which reduces the measured coherence and required the 
calculation of a coherence threshold. This external extrane- 
ous low frequency noise source is the major contributor to 
the three far-held microphone three signal coherence 
shown in Fig. 8. 

Figure 8 clearly shows the broadband noise spectrum 
and the tonal component where it dominates. Above this is a 
broadband spectrum interrupted in places by tonal noise. 
The 1BPF, 2BPF, and 3BPF tones are clearly observable 
along with many other tones in the three signal aligned and 
un-aligned power spectra. Note again, that one can view 
tones approximately 10 dB below the measured auto-spec- 
trum in the three signal coherent aligned and un-aligned 
power spectrum. Similar behavior appears at other turbofan 
engine speeds. 

C. Coherent output power 

Two typical coherent output power 11 ' calculations using 
combustor pressure sensor 1 and 2 with the 120° far-held micro- 
phone (signal source 5) are presented in Fig. 9 for a test condi- 
tion of 1622 rpm (N1 Corr.). These results are shown so that 
they can be compared to auto-spectrum, G55, shown as the top 
curve composed of squares in Fig. 9. To obtain the aligned 
coherent output power using combustor pressure sensor 2 and 
signal source 5 shown by the open diamonds, a time delay of 
0.116875 s was used. To obtain the aligned coherent output 
power for combustor pressure sensor 1 and signal source 5 
shown by the open circles, a time delay of 0.128083 s was used. 
The plane wave combustion noise is apparent between 0 and 



FREQUENCY, Hz 

~ °~" Y 2 ,s I ,=6148/48000=0. 128083 G ss aligned 
I , 6181/48000 0.128771 G JS aligned 
* — y 2 „ I ,=0 G ss unaligned 
""®"”y 2 2s I i=o G SJ unaligned 
G 55 y 2 nn (nd = 250) 

G 5J y 2 i,„(nd-500) 


FIG. 9. (Color online) Total sound power, aligned and un-aligned coherent 
output power calculation for a test condition of 1622 rpm (N1 CORR), using 
signal 1 (pressure sensor 1 at 127° clockwise from top dead center viewed 
from rear) and 2 (pressure sensor 2 at 337° clockwise from top dead center 
viewed from rear) with signal 5 (150 foot microphone at 120°). The light 
dashed line is based on an estimate of the upper 95% confidence limit of the 
un-aligned coherence threshold using nj = 250(yj[ (f, Pj = 0.95, 
n r i = 250) rs 0.012). The heavy line is based on an estimate of the upper 
95% confidence limit of the un-aligned coherence threshold using 
n s = 500(y^(f, P, = 0.95, n s = 500) sa 0.006). 



— Aligned coherence 

Unaligned coherence 

Threshold coherence, 0.012 

Threshold coherence, 0.006 

upper confidence interval 

lower confidence interval 


FIG. 10. (Color online) Aligned coherence, y^y(f), and un-aligned coher- 
ence, |yJ nJn (/)| ? calculation for a test condition at 1622 rpm (N1 CORR) 
using signal 1 (pressure sensor 1 at 127° clockwise from top dead center 
viewed from rear) and 5 (150 foot microphone at 120°). The broken horizon- 
tal line at y J (f, Pi = 0.95, n, j = 250) ^ 0.012 represents an estimate of 
the upper 95% confidence limit of the un-aligned coherence for «^=250. 
The unbroken horizontal line at y\ yn (f, Pi = 0.95, n s = 500) ~ 0.006 rep- 
resents an estimate of the upper 95% confidence threshold limit of the un- 
aligned coherence for n s =500. 


200 Hz in both aligned coherent output power spectra 
yisiflGsstf) | aligned and fls (/) G 55 (/) | aligned ■ Both spectra are 
clearly above their corresponding un-aligned coherent output 
power spectra calculated as yisjunaiignedCss)/) and 
725 1 unaligned C^55 (/") . The un-aligned coherent output power 
spectra were obtained using the as measured signals with no 
time delay. 

Statistical threshold coherent output power curves given by 
yl„y„( n d = 250)G 55 (f) and y 2 Xn y a (n s = 250) G 55 (f) are also 
shown in Fig. 9. The measured aligned coherent output power 
spectra near 100 Flz ( 7,5 (/) G55 (/') | aligned an d yj 5 (f)G 55 
if) (aligned) are clearly above the statistical coherence threshold 
output power curves. These results suggest the presence of indi- 
rect combustion noise at low frequencies due to interaction of 
hot spots with the turbine since they are similar to the results 
discussed by Miles 27-29 for a smaller turbofan engine. 

The results also suggest that combustion noise due to 
higher combustion duct modes is present near 300, 400, and 
600 Hz. Note that the noise floor is about 20 dB below the 
measured auto-spectrum. 

The corresponding aligned and un-aligned coherence 
between pressure sensor 1 and the 120° far-held microphone 
(sensor 5) are shown in Fig. 10. The upper and lower coher- 
ence confidence interval bands for the aligned coherence are 
also shown. In addition, the lower horizontal statistical coher- 
ence threshold line 7 ?„y„( w s) = (/, F/ = 0.95, n s = 500) 

« 0.006 and the higher horizontal statistical coherence 
threshold line 7 2 n>Vi («d) = 7 2 nV;i if- Pi = 0.95, « s = 250) « 0.0 1 2 
calculated from Eq. (7) are shown. The plane wave combus- 
tion noise is apparent between 0 and 200 Hz in the aligned 
coherence function. These are the coherence functions used 
to create Fig. 9. 
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TABLE II. Time series estimation parameters. [r = 48 000 samples/s, NP — number of points in FFT = 4096, 
P — 0.95, B e = r/NP = 11.71 875Hz, n d — B e T tota \ } = number of disjoint (independent) segments] 


Total record 
length, r tota |’s 

Overlap 

n 

V 2 = 

'x n y„ 

i - (i — / > /) 1/( ”— *) 

3-signal noise 
floor, dB 

Coherent output 
power noise floor, dB 

20 

0 

234 

0.0127749 

-9.468 

-18.936 

20 

50% 

S3 468 

0.00639431 

-10.971 

-21.9421 

120 

0 

1406 

0.00221992 

-13.3582 

-26.7164 

120 

50% 

» 2812 

0.00106515 

-14.8629 

-28.7259 

300 

0 

3515 

0.00085215 

-15.3474 

-30.6948 

300 

50% 

« 7030 

0.000426105 

-16.8524 

-33.7048 


XI. DISCUSSION 

Table II shows that using longer record lengths, 7 tot .,|, of 
2 min or 5 min provide a significantly lower noise floor. 
Note that the coherent power noise floor is always less than 
the coherence noise floor based on the three signal method 
being reduced by the square of the coherence. 

Two standard statistical tests were applied to the Fisher 
z-transformed un-aligned coherence results obtained analyz- 
ing simulated time histories and measured turbofan engine 
time histories to determine if they had a normal distribution. 
A norm Q-Q comparison plot showed the simulated time his- 
tory data that produced a more normal distribution than the 
turbofan engine data which had outliers, i.e., points that are 
numerically distant from the rest of the data. The second test 
was to use a histogram plot. This also showed the presence 
of outliers in the measurements made using the turbofan 
engine data. The threshold statistical coherence based on the 
number of independent samples provided a conservative 
value. This value appeared appropriate for the turbofan 
engine data. However, it appeared that a smaller value based 
on the actual number of overlapped samples would be more 
appropriate for the simulated data. 

Finding hidden periodicities in signals with additive 
noise is a signal processing problem with many solutions. 11 
However, the aligned and un-aligned coherence method 
might be useful in cases where the problem involves long 
time histories having an unknown number of multiple peri- 
odic signals and a bad signal to noise ratio. The un-aligned 
coherence procedure has very low computational cost com- 
pared to the other signal processing methods. The method 
might be useful in performing blind test comparisons to vali- 
date and compare these other methods performance on non 
synthetic or simulated time histories. Looking at the literature 
such a method does not appear available. 

The three signal coherence technique calculations 
using one far-field microphone ant two combustor pressure 
sensors shown in Fig. 7 and the total sound power calcula- 
tions shown in Fig. 9 suggest the presence of indirect com- 
bustion noise at low frequencies due to the interaction of hot 
spots with the turbine and are similar to the results discussed 
by Miles 27-29 obtained using a smaller Honeywell TECH977 
turbofan engine. The results obtained in the study of the 
smaller turbofan engine were identified using the convection 
time delay associated with the flow velocity of the entropy in 
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the combustor. For this smaller engine, coherence values 
between the combustor sensor and the far-field microphones 
were much larger and the cross-spectrum phase angle pro- 
vided information on the convection time delay which 
allowed the identification of indirect combustion noise. 

XII. CONCLUDING REMARKS 

In calculating autospectra and cross-spectra using the aver- 
aged periodogram method, the records used are frequently 
overlapped to reduce the variance. In calculating the coherence 
threshold using a statistical method one may use either the 
number of independent samples or a larger number correspond- 
ing to the number of overlapped records used to create the aver- 
age. This paper indicates how applying the Fisher z-transform 
to the un-aligned coherence can aid in making the proper selec- 
tion of samples and produce a reasonable coherence threshold. 
The paper shows that using the number of independent samples 
produces a conservative statistical coherence threshold useful 
for the turbofan engine data studied herein. 

For the simulated time histories based on Gaussian sta- 
tistics used herein, the statistical coherence threshold based 
on a larger number of segments corresponding to the number 
of overlapped records had a lower value. This value is shown 
to be more appropriate by use of the aligned and un-aligned 
coherence and the Fisher z-transform. 

A method has been developed to help identify the coher- 
ent combustion noise source using far-field microphone data 
contaminated with low frequency noise from other sources. 
The method uses an aligned and un-aligned coherence tech- 
nique to determine the significance of the coherence and to 
establish a threshold level for the coherence. This method is 
compared with a statistics based method that can be used to 
calculate the coherence threshold. When used with the clas- 
sical three signal coherence power technique and the coher- 
ent power method these methods provides information on 
coherent combustion noise in the presence of an external ex- 
traneous low frequency coherent noise source. In addition, 
the un-aligned coherence procedure provides a method to 
separate tonal information from random process information. 
The procedure at the same time as it calculates a coherence 
threshold also reveals periodicities, spectral lines, and 
undamped sinusoids hidden or concealed by broadband tur- 
bofan engine noise. Comparison of these results with those 
obtained by another turbofan engine combustion noise 
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source detection study using a smaller engine suggest the 
low frequency combustion noise observed is due to an indi- 
rect combustion noise source. It is expected that this method 
will improve turbofan source diagnostic tests. 
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